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We report on our theoretical study of the inverse TMR effect in the spin polarized transport 
through a narrow channel. In the weak tunneling limit, we find the ordinary positive TMR. The 
TMR changes its sign as the transmission probability becomes large close to a unity. Our results 
might be relevant to the magnetic tunnel junction with a pinhole or a quantum point contact. 
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I. INTRODUCTION 

Recently spin-polarized transport has attracted lots of 
attention because of its potential applications to spin- 
electronic devices. Magnetic tunnel junction (MTJ)ji^ 
consisting of two ferromagnetic (FM) electrodes with an 
insulating layer sandwiched in between two, is one exem- 
plary realization of the spin-polarized transport and show 
a large tunneling magnetoresistance (TMR) effecti^ 
Stable MTJs with large TMR are now routinely fabri- 
cated. The tunneling current is modulated by the rel- 
ative orientation of magnetizations in the two FM elec- 
trodes. According to the Julliere's mode£, the TMR 
ratio is given by 

R P 1 - P L P R v 1 

Here Rp and Rap are the tunnel resistance for the par- 
allel and antiparallel alignment of the MTJ, respectively, 
and Pl and Pp are the electron spin polarizations of two 
FM electrodes. Usually more current flows in the par- 
allel alignment of magnetization such that the TMR is 
positive. 

Recently some experimental groups&LSi2iiS observed 
the inverse TMR effect for which more current flows in 
the antiparallel alignment of magnetizations in two FM 
electrodes than in the parallel alignment. The inverse 
TMR in Co/SrTi0 3 /Lao. 7 Sro.3Mn0 3 is ascribed 6 to the 
selective tunneling of d electrons of Co through the in- 
sulating SrTiC>3 barrier. To investigate the role of the 
oxide barrier on the sign of TMR, MTJs with various 
oxide barriers were studied J. Theoretical studiesii*i2ii^ 
elucidated the important role of the surface DOS and 
the electronic wave functions in the barrier in deter- 
mining the nature of tunneling electrons or the sign of 
the TMR. The surface density of states for Co and Fe 
near the Fermi energy is dominated by the d-band, and 
the s(c?)-band in Co has the positive (negative) spin po- 
larization (SP), respectivelyiii*i2ii^ When the Al oxide 
layer is sandwiched between two Co or Fe metals, the s- 
band electrons relatively easily tunnel through the barrier 
while the d electrons are strongly suppressed in the bar- 
rier. The s electrons are mainly responsible for the tun- 
neling leading to the positive SBii*i2*i^. For the barrier 
of SrTiOa, the opposite is true, leading to the negative 
TMR or inverse TMR. 



The inverse TMR was also observed^ when the spin- 
polarized electrons tunnel resonantly through the impu- 
rity states in the barrier. More current in the antiparallel 
alignment can flow resonantly through a localized state in 
the barrier when the impurity is located asymmetrically 
inside the barrier with respect to two FM electrodes. In 
this case the inverse TMR is determined mainly by the 
structure of the impurity states in the barrier. 

The inverse TMR can occur due to the Kondo effect 
when the impurity spin is sandwiched in between two 
FM electrodes^ The Kondo resonance peak can develop 
even when the magnetic impurity is coupled to the spin- 
polarized electrons of ferromagnetic metals. In the an- 
tiparallel alignment of two FMs, the impurity states of 
spin up and spin down are more symmetrically coupled 
to the corresponding spin states of conduction electrons 
than in the parallel alignment. When the coupling con- 
stants of both spin states are more or less equal, the 
whole Kondo resonance peak develops at the impurity 
siteii - Either in the case of asymmetrical couplings be- 
tween two FMs for the antiparallel configuration or in 
the parallel configuration, the Kondo spin states are split 
by the exchange field of two FMs such that the Kondo 
resonance peak is split into twoii2ii4 Since the Kondo 
resonant state acts as the electric current channel, more 
current can flow in the antiparallel alignment than in the 
parallel alignment, resulting in the negative value of the 
TMR. The inverse TMR was indeed observed in the Ni- 
Cgo-Ni systemiS in the Kondo regime. 

In this paper we propose theoretically another possi- 
ble mechanism for the inverse TMR, based on the simple 
model study. The model system is schematically dis- 
played in Fig. ^ where the spin polarized electrons of 
one FM electrode pass through a narrow channel into 
the second FM electrode. When the transmission prob- 




FIG. 1: Schematic display of nanoscale magnetic tunnel junc- 
tion 
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ability is small between two FM electrodes, the TMR is 
normal or positive. On the other hand the sign of TMR 
can become opposite to the case of weak tunneling when 
the transmission probability is large close to a unity. 

The rest of this paper is organized as follows. In Sec. 
II, the model Hamiltonian is introduced for the nanoscale 
magnetic tunnel junctions and the spin-polarized current 
is formulated using the Keldysh nonequilibrium Green's 
function technique. The results of our work are presented 
in Sec. Ill and a conclusion is included in Sec. IV. In 
Appendix, we present the tight-binding Hamiltonian ap- 
proach to the inverse TMR effect in the nanoscale mag- 
netic tunnel junctions. 



II. FORMALISM 

The considered model system is the magnetic tunnel 
junction with a very small junction area (linear dimen- 
sion is of the order of a few Fermi wavelengths). When 
the junction area is very large compared to the Fermi 
wave length scale, the Slonczewski models will be rele- 
vant. The ferromagnetic (FM) metals, which we conve- 
niently call the left (p = L) and right (p = R) electrodes, 
are described by the two conduction bands of majority 
and minority spins. The relative direction of magneti- 
zations in two electrodes is chosen to be arbitrary. It is 
convenient to write the spin-polarized conduction band 
in the diagonal basis or the spin quantization direction is 
chosen to be the direction of magnetization. 



Hp — 
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ka 
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Here c* pka and c p ka are the creation and annihilation op- 
erators, respectively, for electrons of wave number k in 
the electrode p = L,R, with the spin direction a = ±. 
a = +(— ) means the spin is aligned parallel (antiparallel) 
to the direction of magnetization. e p k a is the energy dis- 
persion relation for electrons in the ferromagnetic metals. 

Since the magnetization directions in the left and right 
electrodes or the spin quantization axes are different, we 
must be careful when writing the coupling Hamiltonian 
between two FM electrodes. The flow of electrons from 
one FM lead to the other is modeled by the transfer 
Hamiltonian liSiil 



H, = 



E E ^Lha,Rk'0c\ kol C R yi3 + H.C. 



(3) 



kk' ap 



The spin quantization axis (the positive z axis) in this 
transfer Hamiltonian is chosen to be the direction of the 
current flow from left to right. In Eq. the spin indices 
a i P —Til- a =T (I) means the spin is directed parallel 
(antiparallel) to the positive z direction. Note that we 
are completely free to choose any direction as the spin 
quantization direction for this transfer Hamiltonian. Af- 
ter fixing the positive z axis, we can select the orthogonal 



x and y axes, perpendicular to the z axis, at our dis- 
posal. Once the coordinate system is fixed, the direction 
of magnetization in each FM electrode can be specified 
by two angles (9 p ,<p p ) in the spherical coordinates with 
p = L,R (left, right electrode). Here 9 p measures the an- 
gle between the magnetization direction and our chosen 
z axis, and (j) p is the azimuthal angle of the magnetiza- 
tion measured from the x axis. thku^Rk'p in Eq. 10 is 
the hopping integral from the spin band f3 in the right 
electrode to the spin band a in the left. At this point we 
consider the most general form of the hopping integral 
which allows the spin flips. The first term in Eq. @ rep- 
resents the transfer of electrons from the right electrode 
(R) with the spin direction (3 to the left (L) lead with 
the spin direction a. This transfer Hamiltonian is rele- 
vant when the transfer of electrons from one lead to the 
other occurs dominantly through a very narrow channel 
(one transport channel). 

It is straightforward to show that the electron annihi- 
lation operators in two different bases of Eqs. © and © 
are related to each other by the equation 
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Let us introduce the electron spinors ($ and \&) and the 
unitary transformation matrix (U). 
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p^r sin 9 -f cos % 
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The electrode Hamiltonian, Eq @, is diagonal in the $ 
basis, but not diagonal in the \& basis. The two bases 
are related by the equation, ^ p k = U p § p k, which is none 
other than Eq. (0J. Using the spinor operators, the model 
Hamiltonian can be written as 



H p 
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(6b) 
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Here we introduced the hopping matrix VL k ,Rk' ■ 

We are now in a position to derive the expression of 
spin-polarized current for our model system. Using the 
nonequilibrium Green's function methodr^^the electric 
current flowing from left to right can be computed by 
thermally averaging the current operator 20,21 



/ = 



P c Lka C Rk>f3 



H.c. 
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It is now more convenient to write the current operator 
in the matrix form. 



I = 



kk' 



(7b) 



Let us now introduce the time-ordered mixed Green's 
functionsi2*i^ which are defined by 

ihG L R(kt,k't') = (TV Lk {t)* Rk ,(t')), (8a) 
ihG RL {kt,k't') = (T* Rk (t)tf Lk , (*')), (8b) 

where the symbol (A) means the thermal average of 
A and T inside the thermal average means the time- 
ordering in the Keldysh contour. In terms of the above 
mixed Green's functions, the current can be written as 

/ = 2eJ2 l ^[VLk,Rk'G< L (k't,kt)] 

kk' 

= f [ delimit [V LktRk ,G< L (k',k;e)] . (9) 

J kk' 

The second line is obtained from the first line after the 
Fourier transform from time to energy variables. 

Since we are mainly interested in the linear response 
conductance of the system, the physical properties are 
determined by the energy structure near the Fermi level 
at low temperature compared with the Fermi tempera- 
ture. This is the case we are going to consider below. 



Since the transport properties are mainly determined by 
electrons near the Fermi level, the dependence of the hop- 
ping integral on the wave vector may be neglected and 
the hopping integral can be replaced by its value at the 
Fermi energy. From now on, the dependence of the hop- 
ping matrix Vlr on the wave vector will be suppressed 
under this approximation scheme. In the tight-binding 
Hamiltonian approach, the hopping matrix Vlr does not 
depend on the wave vector (see the Appendix for details.) 
Under this approximation, it is convenient to introduce 
the wave- vector-summed Green's functions^ to simplify 
the algebra. 



G P (M') = ^G p (fc;M'), P = L,R, (10a) 



ihG p (k; t, t') 
Grl (t, t') 



(T* pk (t)*l k (t')), 
J2G RL (kt;k't'). 



(10b) 
(10c) 
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Note that G p (k; t, t') is the Green's function of electrons 
in each lead p — L, R when the coupling between two 
FM electrodes is absent, while GRL(kt; k't') is defined in 
Eq. I)8bj) . Using the Feynman diagrams, we can readily 
show that the mixed Green's function is determined by 
the Dyson-like equation^ 



G RL (t,t') = G RL0 (t,t') + f dhGRLofot^VLRGRUtutf), (11) 

Jc 

where the auxiliary Green's function Glro is defined as GRLoit.t 1 ) = J c dtiGR(t,ti)VRLGL{ti 7 t'). Here C in the 
integral sign is the Keldysh contour embracing the real time axis. The Green's functions Gl /r in the above equations 
can be interpreted as the Green's function of an electron located at the local tip site of the left/right lead in the 
absence of the coupling between two FM leads (see the Appendix.) After analytically continuing— to the real time 
axis from the Keldysh contour, it is straightforward to find the lesser Green's function G RL 

G< L (e) = [l-G^ (e)U Lfl ]' 1 G< L0 (e)[l-U L flG?j L0 (e)]- 1 , (12) 

where the axillary Green's functions are 

G r R a L ,{e) - G^(e)VR L Gr(e), (13a) 
Gi L0 (e) = G<(e)V RL G a L (e) + G R (e)V RL G<(e). (13b) 

The G <,r,a denote the lesser, retarded, and advanced Green's functions, respectively. 

To find the expression of the spin-polarized current, we need the explicit forms of the lesser, retarded, and advanced 
Green's functions for electrons in the FM electrode. We have to note that the Green's functions of the FM electrodes 
in the above equations are not diagonal, because we worked in the spin basis of the transfer Hamiltonian. 

ihG p (k;t,t') = (T* pk (t)*l k (t')) = U p (T$ pk (t)$l k (t'))U;, (14a) 
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Here g p is the Green's functions of the FM electrodes in the spin diagonal basis and related to G p by the unitary 
transformation, G p — U p g p U^. Using the above relations, the desired Green's functions can be found for the model 
Hamiltonian of Eq. J2J 

OS-W = U p J2^-^ ±iS ^ [e _ ep ° ±i ^)u}, (15,) 
G<( £ ) = 2,f r W t J2 ( S(e {(e _° epk _)) Oj- (15b) 



Here an infinitesimally small positive number S is in- 
cluded for the definition of the retarded (+ sign) and ad- 
vanced (— sign) Green's functions, G p ' a . f p (e) = f(e—fx p ) 
is the Fermi-Dirac thermal distribution function in the 
lead p — L. R and \i p is the chemical potential shift of 
each lead caused by the source-drain bias voltage. Let us 
define R p ± and N p ± as 



R p± (e) 
N p ±{e) = ^£(e- e pk ±), 



(16a) 
(16b) 



where V means the principal part of the integral. N p ±(e) 
is none other than the density of states (DOS) for 
the majority/minority spin in the FM lead p = L, R. 
The retarded, advanced, and lesser Green's functions 
[Eqs. (|15a|) and l|15b(l ] can be written as 

G r p a (e) = U p [RpTiTrN p ]Ul, (17a) 
G<(e) = 2*f p (e)U p N p U}, (17b) 



where R p = 



(Yv) u*lNp= ( Np o + nI)- When 
the Fermi level lies inside the band, R p ±(e = Ep) is ex- 
pected to be small due to the mutual cancellation of elec- 
tron (above the Fermi level) and hole (below the Fermi 
level) contributions. In this case or in a wide conduction 
band limit^ we can set R p ±{e = Ep) ~ 0. The retarded 
and advanced Green's functions are given in simple forms 
in this case. 



Gr,a 
p 



(e) 



TinU p N p Ul. 



(18) 



For half metals with only one spin band occupied at the 
Fermi level, the value of R for unoccupied spin band may 
not be neglected such that the retarded and advanced 
Green's functions may have real component. 

When the Fermi level lies well inside the band for both 
spin directions or in a wide band limit^ we can approxi- 
mately set R p ± as 0. In this case, inserting the Eqs. 117bll 
and l|18fl into Eqs. i|13|) . the expressions of G 
Grlo can ^ e simplified as 



RLO and 



-tt 2 N r V rl N l , (19a) 
2t7T 2 [f R (e) - f L (e)] N R V RL N L , (19b) 



where N p — U p N p U^. Inserting these two Eqs. (|19|l into 
Eq. Q12[L we find for the lesser Green's function G^ L 



G< L = 2tn 2 [f R (e)~f L (e)} 



XNrVrlNl 



1 + ^NrVrlNlVlr 

1 



1 + k'VlrNrVrlNl 



(20) 



Combining this lesser Green's function and Eq. @, the 
expression of the spin-polarized current can be written 
as 



I = — 



delmTr [V LR G< L (e)] 
de[f R (e) - f L (e)]T(e) 



where 



T(e) = 2ReTr[l + r]" 1 r[l + r]" 

r = tt 2 NrVr L n l v L r. 



(21) 

(22a) 
(22b) 



In obtaining the above form of T(e), the trace property of 
the matrix product was invoked. Here T(e) is the trans- 
mission probability through a narrow channel for spin- 
polarized electrons, when the magnetization directions of 
two FM electrodes are arbitrary. The information about 
the orientation of two magnetizations is included in the 
DOS matrices, Nl and Nr. The matrix T measures the 
effective transfer rate of spin-polarized electrons from one 
FM electrode to the other FM electrode, and contains the 
dependence on the relative orientation of two magnetiza- 
tions. The above form of the transmission probability is 
valid for the wide band case or when the Fermi level lies 
well inside the bands for both spin directions. The linear 
response conductance G is then given by the expression, 
G = =f- ■ T(e — 0). Below we are going to confine our 
interest to the linear response regime and to the case of 
a wide band limit. 

We can give some physical meaning to the density ma- 
trix N p — UpNpUp'. After some algebra, N p can be ex- 
pressed as 



N„ = 



N, 



N p0 [l + P p fip ■ a], 
1, 



p0 



Pn 



-[N p+ +N p _ 



N 



Np^ 



N n 



(23a) 
(23b) 

(23c) 
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Here a is the Pauli matrices, N p o is the average DOS 
of the majority and minority spin bands. P p measures 
the spin polarization at a given energy e of the FM elec- 
trode p = L, R and agrees with the definition of the spin 
polarization in the Jullicrc model. h p is the unit vec- 
tor defined by the two angles, 8 P and <p p in the spherical 
coordinate system and is aligned parallel to the magne- 
tization direction of the p = L,R FM electrode. 



tLB.8 a f3- Electrons do not flip their spin while passing 
through the narrow channel. In this case the expression 
of r in the transmission probability, Eq. I|22a|l . is simpli- 
fied as 



r = 7r 2 \t LR \ 2 N R N L 

= 7 (1 + P R h R ■ a)(l + P L h L ■ a). 



(24) 



III. RESULTS AND DISCUSSION 

In this section we study the TMR behavior for the 
wide band case with the transfer Hamiltonian: t^^ jip — 

1 



Here 7 = ir 2 NLoN R o\tLR\ 2 measures the average transfer 
rate of electrons for two spin directions. The expression 
of the transmission probability T is 



T(e) = 4 7 x 



1 + Pl ■ Pr + 7(1 ~ ~ Pr)(2 + 7 + lP L ■ Pr) 

[1 + 2 7 (1 + Pl • Pr) + 7 2 (1 - P 2 L )(l - P 2 R )] 2 

I 



(25) 



Note that Pl = ulPl and P R = n R P R . This is the 
most general result for the magnetic tunnel junctions 
with a narrow junction area. Note that the conductance, 
2e 2 /h x T(0), depends only on the relative angle of two 
magnetization directions. In the weak tunneling limit 
or 7 -C 1, the transmission probability can be approx- 
imated as T = 4 7 (1 + P L ■ Pr) = 4 7 (1 + P L P R cos9). 
The conductance is then given by the expression, G = 
G (l + P l Prcos9) with G = 8je 2 /h. This is the fa- 
miliar expression of the conductance which is in agree- 
ment with the Jullicrc model. 5 The TMR ratio is then 



TMR = 2P L P R /(1-P L P R ). In the weak tunneling limit, 
the TMR is determined solely by the polarization of the 
ferromagnetic metals. In the strong tunneling limit the 
angle dependence of the conductance is still an even func- 
tion of 6, but is more complicated as shown in Eq. 125|l . 



When two 
(N L+ = N R+ 



ferromagnetic electrodes are 



N + and N L - 



N 11- 



th.e 



same 
with 



the possible difference in the magnetization direction, we 
have Pl = P R = P- The expression of T, Eq. l(25fr. is 
further simplified 



J 



T(cos( 



= 4 x 



7(1 + P 2 cos6>) + 7 2 (I - P 2 ) 2 {2 + 7 [1 + P 2 cosl?]) 



[1 + 27(1 



P 2 cos ( 

I 



+ 7 2 (1-P 2 ) 2 ]' 



(26) 



Especially for the parallel (cos 9 = 1) and antiparallcl 
(cos 6 = — 1) configurations, the transmission probabili- 
ties, respectively, given by the expressions 



Tp 
Tap 



274 



27- 



(1 + 7+) 2 (1 + 7-) 2 ' 
V7+7- 



(27a) 
(27b) 



Here 7± is defined by the equation r y± = Tr 2 N±\tLR\ 2 (0 < 
7± < 1), and is related to the transmission probability 
T± by the relation 



4 7± 



(l+7±) 2 



(28) 



T+ is the transmission probability for electrons from the 
majority band in the left electrode to the majority band 
in the right electrode, while T_ from minority to minor- 
ity. Note that the average transfer rate 7 and the po- 
larization P are related to the transfer rates of majority 
and minority spins by the equations 



7 
P 



(29a) 
(29b) 



The transmission probabilities, Tp 
expressed in terms of T± as Tp = 
T A p = 4^7^(1 + ^T+)(l + ^l~TZ) [j%TZ 



and Tap , can be 
i(T+ + T_) and 
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FIG. 2: Inverse TMR. (a) Dependence of TMR on the trans- 
mission probabilities of electrons for majority (T+) and minor- 
ity spin (T_). The dashed line is the boundary between the 
positive and negative TMR. (b) Sign of TMR in parameter 
space of (T+ , T_ ) . TMR vanishes along the solid and dashed 
lines. Two lines meet at T+ = T_ = 2/3. The TMR becomes 
negative when the tunneling probabilities for electrons of both 
spin directions is close to a unity. The curve T+ = T_ (solid 
line) corresponds to the nonmagnetic electrodes. 



(l + yi-rVxi + yi - iqr 2 . 

In the weak tunneling limit T± <C 1 or "f± <C 1, the 
Eq. I|2()l) can be approximated as T(cos6) — 47(1 + 
P 2 cos 9) . Especially for the parallel and antiparallcl 
configurations, we have Tp = 2(7+ + 7_) and Tap — 
4^/74.7- . Since 7+ + 7- > 2^/7 + 7_, we have the con- 
ductance inequality Gp > Gap- This inequality rela- 
tion is true generally when the transmission probability 
is weak. That is, the current flows more in the parallel 
configuration than in the antiparallel configuration. 

In the strong tunneling limit, the above inequality in 
conductance can be reversed and more current can flow in 



FIG. 3: Dependence of conductance on the relative angle be- 
tween two magnetization directions with varying transmission 
probability T+ for the majority spin. In panel (a) P = 33% is 
chosen to simulate the Ni electrodes;— while P = 55% in panel 
(b) for CoFe electrodes. The solid line (T+ << 1) corresponds 
to the Julliere result or the very weak tunneling limit. 



the antiparallel alignment than in the parallel alignment. 
To see when the current flows more in the AP than in 
the P, we study T P and T A p [Eqs. (|27ajl and (|27E|l ] in 
the plane of (7+ , 7_ ) or (T+ , TL) . The curves on which 
Tp = Tap are given by the equations 



= 7+ - 7- 



(30a) 



= ( 7+ 7--l) 2 -2V7^(l+7+)(l+7-)(30b) 

Along these two curves the conductance in the P and 
AP configurations is equal. The curve 7+ = 7_ corre- 
sponds to the case that Tp = Tap and the magnetic 
polarization at the Fermi energy is zero or P = 0. The 
second curve, Eq. I|30b|l . is the boundary between two 
regions: Tp > Tap (close to the origin) and Tp < Tap 
(strong tunneling regime), and corresponds to the case 
that P ^ 0, but Tp = Tap- Two curves meet each other 
at (7+, 7-) = (2 — v3, 2 — V3), which corresponds to 
the point (T+, T_) = (2/3, 2/3) in the plane of transmis- 
sion probabilities. The behavior of TMR is summarized 
in Fig. 12(a) and (b). The value of TMR is plotted as 
a function of (T+ , T_ ) in the panel (a) and the bound- 
ary between normal and inverse TMR is marked by the 
dashed line. The panel (b) provides the TMR diagram 
in the (T + ,T_) plane. Especially in the strong tunneling 
limit, the inverse TMR can be realized in the nanoscale 
magnetic tunnel junctions. 

Fig. |3| displays the conductance ratio, G(8)/Gp, as a 
function of the relative angle between two directions of 
magnetization [Gp — G(8 = 0)]. While the polarization 
is fixed, the transmission probabilities are varied. In the 
weak tunneling limit the solid curve traces the cosine 
function. With increasing transmission probabilities, the 
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curves deviate from the Julliere result and furthermore 
the conductance minimum can change into a maximum 
leading to the inverse TMR. 

IV. SUMMARY AND CONCLUSION 

To summarize, we studied the spin-polarized transport 
through a narrow channel using the transfer Hamilto- 
nian approach and the nonequilibrium Green's function 
method. Our analysis applies to the case when the Fermi 
level lies well inside both majority and minority spin 
bands or to the wide band limit. In the weak tunnel- 
ing limit we reproduce the Julliere results and the TMR 
is positive. On the other hand, in the strong tunneling 
limit the inverse TMR opposite to the weak tunneling 



case is found. We presented another theoretical scenario 
for the inverse TMR effect. Our model study may be 
relevant to the magnetic junctions with a quantum point 
contact between two ferromagnetic metals or to the thick 
planar MT Js with a pinhole which accommodate only one 
transport channel. 
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APPENDIX A: TIGHT-BINDING MODEL FOR THE FERROMAGNETIC ELECTRODES 



To illustrate the relevance of our approach to real systems, we consider the description of the ferromagnetic (FM) 
electrodes in terms of the tight-binding Hamiltonian. Two FM electrodes are coupled by the hopping term between 
two atomic sites. Though our formalism can be generally applied to any types of lattice, we confine our study to the 
lattices for which the analytic solution of the Green's functions can be readily found. The examples include the Bethe 
lattice out of which the simplest one is the semi-infinite one-dimensional tight-binding atomic chain (see Fig. El) 

The left and right electrodes are described by the tight-binding Hamiltonian 

Hl = E E e ^ c iW -EE [W L *4 +Ia c la + H.c] , (Ala) 

i<-lcr=± i<~2<J=± 

^ = EE ^CjW "EE [ W R-4+la^a + H. C ] . (Alb) 
i>l cr=± i>l cr=± 

The left (right) FM electrode is terminated at the site of i = —1 or i = L (i = 1 or i = R), respectively. The 
i = L site will be called the left terminal site and the i = R site the right terminal site. The electron creation (cj^) 
and annihilation (c^) operators are defined in the spin diagonal basis, a = +(— ) means the spin is aligned parallel 
(antiparallel) to the magnetization, respectively. However the spin quantization directions of two leads are different 
due to different magnetization directions. The on-site energies, tLa and eji a , include the exchange splitting between 
two spin directions as well as a possible voltage difference between two FM leads. In order to take into account the 
different bandwidth of spin-up and spin-down bands, the hopping integrals, 24 Wlct and Wn a , include the spin index 
a = ±. Magnetic polarization is determined by the values of on-site energies and hopping integrals. 

Adopting the same convention for the direction of spin quantization (chosen to be the direction of the current flow 
from left to right) as in the main text, the hopping term which is responsible for the flow of electrons between two 
leads can be written as 

Hl = E [tLcRffciaCRp + ff.c] . (A2) 

a/3 

Here a., (3 =],[ and the ] {{) direction is parallel (antiparallel) to the spin quantization direction. Since the spin 
quantization directions in Eqs. (|A1|) and Eq. (|A2(I are different, the electron operators in two spin bases are related 

w La t W Ra 

• — • — • — • o — o — o— 

-4 -3 -2 -1 (L) 1 (R) 2 3 



FIG. 4: Tight-binding model for the FM electrodes. The FM electrodes are described by the semi-infinite chain. The site of 
i = — l(i = 1) is the terminal site of the left (right) electrode, and is denoted as i = L(i — R) for the notational convenience. 
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to each other by the unitary transformation. For details, see the main text. 

(A3a) 
(A3b) 

The two bases are related by the equation, \ff p = U p $ p for p = L,R. Obviously the two bases at other sites are 
related to each other by the same unitary transformation. Using the nonequilibrium Green's function method, the 
spin-polarized current can be expressed as 




■ J de JmTrV LR G< L (e), (A4) 
where the hopping matrix Vlr and the Green's function are defined as 



V LR = [ l Ltm \ LhRi ) , (A5) 



ihG RL (t,t') = (T* R (t)tf L (t')). (A6) 

Note that the hopping matrix has no dependence on the wave vector in the tight-binding Hamiltonian approach. Here 
GnL(t,t') is the time-ordered Green's function and G^ L is the lesser part of G R l- After Fourier transform over time 
variable, the lesser Green's function G^ L can be found in a closed form (see the main text and also |2(j.) 

Gi L (e) = [1 - G^WlrV 1 G< L0 (e) 

x[l~V LR G RL0 (e)}-\ (A7) 

GTlo(£) = G r A a (e)V RL G2 a (e), (AS) 

Gi L0 (e) = G<(e)V RL Gl(e) + G R (e)V RL G<(e). (A9) 

Here Gl and G R are the Green's functions at the left and right terminal sites, respectively, when the coupling between 
two FM electrodes is absent. 

ihG p (t,t') = (T* p (t)*t(f)) = U p {T%{t)^ p {t'))Ul p = L,R. (A10) 

The superscripts r, a and < mean the retarded, advanced, and lesser part of the Green's function. To compute the 
spin-polarized current, we have only to find the Green's functions at the terminal site (L or R) when the coupling 
between two FM electrodes is absent. It is convenient to define the terminal Green's function in the spin diagonal 
basis 

«*fcO - (r* l(i) *»(0) - (<^+«>4 + «'» {Tcr J 4 _ {l , }) ) ■ (An) 

This Green's function g p = (^ 9p + g ° ^ is diagonal and is related to G p by the unitary transformation, G p = U p g p U^. 

At this point we had better take a concrete example for the terminal Green's functions, and g R (in the spin 
diagonal basis). As an example we consider the semi- infinite chain model for the FM electrodes. It is relatively simple 
to find the self-energy of the terminal Green's function, when we take into account the repeating structure of the 
chain. The self-energy can be expressed in terms of the Green's function itself (p = L.R). 

v a +>) _ (\W p+ \ 2 g p+ {t,f) \ , . 

MM) " 1. \W P A 2 9 P -{t,t'))- {AlZ) 

From the Dyson equation, we can find the desired Green's functions in a straightforward manner. The retarded and 
advanced Green's functions are 

g v°^ = a. ■ IaZ ia t =V- (A13) 

e - e pCT ± iy/4\W pa \ 2 - (e - e pCT ) 2 
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The retarded part (r) corresponds to the + sign, while the advanced part (a) to the — sign. Quite in general, the 
retarded and advanced Green's functions at the terminal site can be written in a form g pa {t) = R pa {e) — iTrN p(T (e) 
and (7p CT = [g p<7 ]* , respectively, such that we can write the matrix Green's function as 



f p a 



(e) = [\® ? (£) iA)) = M*)Ti«N p (e). 



Here N pa is the local density of states (DOS) or the local spectral function at the terminal site in the spin state 
a = ±. The lesser Green's function at the terminal site is simply related to the DOS matrix N p — ( q + ^ ) by the 



expression, g< = 2-k f p (e)N p . 

Once the terminal Green's functions are known in the diagonal basis, the algebra becomes straightforward to reduce 

the expression of the spin-polarized current, Eq. (|A4|) . Then the desired Green's functions at the terminal sites are 
(p = L,R) 

G r p a (e) = U p ^+ (e) £° (e) ) Ul (A15) 

G<(e) = 2vf p (e)U p N p V$. (A16) 

To simplify the notations we rewrite 

G r p a (e) = R P T^N P , (A17) 

G<(e) = 2nf p (e)N p , (A18) 

where N p = U p N p U^ and R p = U P R P U^,. The expression of G^ L0 can be rearranged as 

G< L0 = 2iir 2 [f R - f L ]N R T RL N L 

+2tt [f R N R T RL R L + /lRrTrlN^ . (A19) 

From now on we consider the case of no spin- flip in tunneling of electrons between two FM electrodes or T R l = t R L^-- 
We also confine our interest to the case that two FM electrodes are the same material. Then we can derive the simpler 
relation for the spin-polarized current. 

G< L0 = 2it RL ir 2 [f R - f L ]N R N L 

+2Trt RL [f R N R R L + f L R R N L ~\ , (A20) 

^{TrlG< l } = t LR Tt{[i-t LR cr RIJ0 ]- 1 Gi L0 

xll-WWr 1 }. (A21) 

From these two relations, we find the simple expression for 1tii£t:Ti jR G < ^ l in the parallel and antiparallel alignments 
of magnetizations such that the expression of the spin-polarized current is 

1 = Tj de V*®-hW T &> (A22) 
T(e) = 2n 2 \t LR \ 2 Tr{[l-t LR G RL0 ]- 1 N R N L 

[1-W^Lor 1 }- (A23) 
Since we are considering two FM electrodes of the same material, we have for the parallel alignment of magnetizations 

Nl = Nr = ^ £) , (A24) 

G RL0 - t R L ^ Q {r _ t . nN _ )2 J . (A25) 

The transmission probability for this case is 

T P (e) = V (A26) 



<7 = ± 
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where g a = + iirN a . For the antiparallel alignment of magnetizations, we have 



N_ 
iV+ 

Gfllo = ^(-R + TwJV+)(i?-TwiV_)l, 



such that the transmission probability is given by the expression 

4ir 2 \t LR \ 2 N+N- 



T A p{e) 



\l-\t LR \ 2 g + g. 



(A27) 
(A28) 

(A29) 



We are now in a position to discuss any possible inverse tunneling magnetoresistance effect for the narrow channel 
magnetic tunnel junctions. 

In a weak tunneling limit or when |tz,.Rg£| <C 1, the transmission probabilities for the parallel and antiparallel 
alignments of magnetizations can be approximated as 



T P = 2it 2 \t LR \ 2 [N* +Nl'_ 



4ir 2 \t LR \ 2 N+N_. 



(A30) 
(A31) 



Obviously Tp > Tap such that more current flows in P than in AP, and the positive TMR ratio is obtained with the 
familiar Julliere expression. 

In a strong tunneling limit, let us see if there is a possibility of Tp < Tap- To be concrete let us consider the 
semi-infinite chain model. Note that the Eq. (|A13|I can be written as 



per 



2W, 



r? 



pa 



4|WV| 2 



(A32) 



and g pa = [gp a ]* ■ Since we are interested in the linear response transport, we have only to consider the terminal 

Green's functions at the Fermi level or e = 0. In a wide band limit or when \e pa \ <C |Wpo-|, the terminal Green's 

functions are reduced to g prT ~ —i/\W pcr \ and <7p CT ~ i/\W pa \ such that our discussion corresponds exactly to our study 

in the main text. When R± = 0, the transmission probabilities for the tight-binding model are exactly the same as 

those in the main text. It is expected that the inverse TMR is possible when \R\ <C N. 
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